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A (globally) neutral two-body system is supposed to obey a pair of coupled Klein-Gordon equations in a con- 
stant homogeneous magnetic field. Considering eigenstates of the pseudomomentum four-vector, we reduce 
these equations to a three-dimensional eigenvalue problem. 

The frame adapted to pseudomomentum has in general a nonvanishing velocity with respect to the frames 
where the field is purely magnetic. This velocity plays a crucial role in the occurance of motional terms; these 
terms are taken into account within a manifestly covariant framework. 

Perturbation theory is available when the mutual interaction doesnot depend on the total energy; a 
weak- field- slow-motion approximation is more specially tractable. 

PACS 03.65.Pm (relativistic wave equations) 



1 



1. INTRODUCTION. 

Prom a relativistic point of view, the constant magnetic field has this peculiarity that it does not 
correspond to a unique laboratory frame. When a constant homogeneous electromagnetic field is seen as 
purely magnetic in some frame (conventionally referred to as laboratory frame), such a frame cannot be 
imique [1] . The directions eligible for the time axis of a lab frame span a 2-dimensional hyperbolic space 
(E^) (longitudinal space). Energy, defined as time component of the momentum, is conserved but it is 
affected by this ambiguity about the lab frame. Therefore a manifestly covariant treatment is specially 
relevant. 

Another important feature of the constant magnetic field is the conservation of pseudomomentum, which 
can be checked in several cases of interest. In these cases, the true momentum (associated with translations 
in configuration space) is not fully conserved, because application of a magnetic field spoils translation 
invariance in the Klein-Gordon operator, but it turns out that the so-called twisted translations [2] operating 
in phase space, preserve the dynamics, which results in the occurence of a conserved vector [2-5] which 
reduces to ordinary momentum in the no- field limit. 

At least this statement is true for n charges in the Galilean theory, supposing that the mutual interaction 
between two charges corresponds to a central force [2]. 

It is also true for one-body and two-body relativistic systems, under rather general assumptions. Check- 
ing the one-body case is straightforward: the four-velocity being p — eA, and the gauge being A= ^q- F, it 
is a mere exercise to verify that C = p + eA commutes with the Klein-Gordon operator. 

For scalar particles, the two-body problem can be posed in terms of a pair of coupled wave equations 
[6,7]. 

In the absence of external field, the dynamics of the system has all the Poincare symmetries, in particular 

it is invariant under ordinary translations; as a result total momentum P is conserved. Then, imposing sharp 
values to its components and after elimination of the relative time, one is left with a reduced wave equation 
involving three degrees of freedom only. At least when the mutual interaction potential does not depend on 
the total energy [8], this reduced equation can be seen as a standard eigenvalue problem (in terms of the 
energy of relative motion). As the constituent masses have been fixed from the outset, solving this equation 
yields a spectrum (generally with a discrete part). Admissible values of are correspondingly selected. 

When an external electromagnetic field is applied to the system, a first problem consists in keeping the 
wave equations consistent with one another. When the field can be seen as purely magnetic and constant (in 
space and time) for some observer, the compatibility requirement, combined with symmetry conditions and 
the demand of reasonable limits when mutual (resp. external) coupling is turned off, is satisfied in closed 
form by an ansatz [1] [9] [10] , provided the that the mutual interaction term is known in the isolated system. 

Now the customary momentum P is not any more a constant of the motion. But we have previously 
proved this result [1,10]: for a large class of mutual interactions (sufficient for all pratical purposes) the total 
pseudomomentum of (globally) neutral systems (that is ei = — e2 = e) is actually conserved. Moreover the 
four components of this vector commute among themselves. This result seems to admit an extension to the 
case of an arbitrary constant electromagnetic field [10] , but only the pure electric or magnetic cases lead to 
a rigorous statement and to wave equations written in closed form, after a suitable transformation [11]. 

It is possible to take advantage of pseudomomentum conservation in order to reduce the degrees of 
freedom, in a way which parallels the customary separation of center-of-mass variables usually performed 
in isolated systems. In addition, a further degree of freedom can be separated out. In the magnetic case 
(considered throughout this paper), this degree of freedom is a timelike variable identified as relative time 
in one of the possible lab frames. 

This situation provides a description of neutral bound states in a magnetic field. Indeed all the com- 
ponents of pseudomomentum can be simultaneously diagonalized. In a generalized sense, the motion of 
the system as a whole is separated out, and inner motion can be considered; in the same spirit as in the 
nonrelativistic theory, it is natural to associate a spectrum to this (pseudo) relative motion. 

Solving the reduced eigenvalue equation will determine admissible values for the square of pseudomo- 
mentum; having bound states in mind, we shall be more specially interested in the discrete spectrum. 

Since pseudomomentum is the generalization of momentum in the presence of field, its square plays the 
role of an effective squared mass and should be observable. We are led to investigate how much this quantity 
is shifted by the field from the total squared mass we would have obtained for the isolated system. 
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In a siiitablc representation, wo shall explicitly write down and reduce the wave equations. 

As the reduced equation is a three-dimensional problem, it is tempting to apply the standard methods 
of perturbation theory. But, even in the simple case where the unperturbed motion is ruled by a mutual 
interaction term which does not depend on P^, a complication arises due to the presence of external field: 
the reduced equation nonlinear ly depends on the total energy. As a result, it is not an eigenvalue equation in 
the usual sense. Fortunately, a spectral theory for energy-dependent perturbations of a standard eigenvalue 
equation has been already presented in the literature [12]. 

Let us now say a few words about the geometric elements corresponding to a constant magnetic field in 
the presence of a constant timclike vector fc, eigenvalue of pscudomomentum. 

A constant magnetic field F'^'' provides a unique and relativistically invariant decomposition [1] of any 
four- vector ^ into longitudinal and transverse parts, ^ = + ^t- 

In particular k = k^ + kT- Assuming that k is timelike, it is clear that kr can accidentally vanish whereas 
fci never does. In fact fci is always timelike. Since fci G (El) it can be considered as defining a preferred 
lab frame, we shall refer to it as the special lab frame. Let us emphasize that this preferred lab frame also 
depends on the state of motion of the system (in contradistinction to the unique lab frame associated with 
an inhomogeneous field). 

It is noteworthy that = and therefore Ik^] is nothing but the energy seen in the special lab frame. 

In general the frame adapted to k (which differs from fc^,) cannot be a lab frame. In view of results 

about pseudomomentum in the Gallilcan theory, it is reasonable to admit that k carries information about 
the motion of the system as a whole. So we call pseudo-rest frame the frame adapted to k. 

\kJ 

With respect to the special lab frame, it has a velocity characterized by the motional parameter e = -rr—r- 

\kL\ 

In the particular case where kj. — 0, then fci, and k coincide and e vanishes; there exists a frame adapted to 
k where F is purely magnetic and the subsequent developments get drastically simplified. Otherwize, terms 
involving the contraction k ■ F arise. We notice that {k ■ F)°' = \k\E°' where i?" is the electric field "seen" 
by an observer moving with momentum k°' {motional electric field). 

This paper is organized as follows. In next section we recall the ansatz which provides explicit equations 
of motion, and simplify a quantity Z which was an essential ingredient in the formulas of rcfs.[l, 10]. In 
Section 3, several expressions given in compact form in that previous work are explicitly developed for 
applications. Section 4 is devoted to the three-dimensional reduction, with emphasis on the role of relative 
energy. Various terms present in the reduced equation are discussed and ordered according to the powers 
of field strenght and motional parameter. In Section 5 we consider normalization and the possibility of a 
reliable perturbation treatment. Last Section is devoted to a few concluding remarks. 

2. BASIC EQUATIONS 

A system of scalar particles can be described by a pair of coupled Klein- Gordon equations = 
■^'ma ^ a,b = 1,2 where ^ has two arguments qi,q2 running in spacetime. We cover all cases of 

practical interest assuming that 

Ha=Ka + V 

In this formula 2Ka is the squared-mass operator for particle a alone in the magnetic field, and F is a suitable 
modification of the term V^^^ which describes the mutual interaction in the absence of external field; more 

generally, the label (0) refers to the no-field limit of any quantity. 

In order to be more specific we separate canonical variables in two classes as follows. 

P = Pl+P2, Q=^(?l+92) 

z = qi-q2, y=^{Pi-P2) 
This provides two sets of standard commutation relations. Useful quantities are 

u' = y'-{yPr/P' 
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The latter is an essential ingredient of mutual interactions; but in order to avoid denominators in calculations, 
it is convenient to employ 

Z = z'^P'^-{z-Pf (2.1) 

We shall assume that 

yW=/(Z,p2,t/.P) (2.2) 
Since our system is globally neutral, it turns out that pseudomomentum is 

C = P+^z-F 

In canonical variables we can separate the transverse pieces from the longitudinal ones, for instance P = 
Pt + Pl, etc. Transverse and longitudinal variables mutually commute. 

Compatibility requires that [Ki — K2,V] =0. Adapting Bijtebier's method [9] to the peculiarities of 
the constant magnetic field, we have performed a canonical transformation say 

= exp{iB) O' = exp{iB) O exp{-iB) (2.4) 

where O is any operator. B was choosen [1,13] such that transformation (2.4) yields 

K[-K^ = y^- P^ (2.5) 

thus compatibility is satisfied through the ansatz [1] 

V' = f{Z,P',y,.P,) (2.6) 

where Z = Z'^^^ = {Z')f=o (it turns out that Z commutes with y^, -Pj,). The explicit form of Z was calculated 
in ref.[l] 

Z ^ Z + 2{Pl z-P-P^ Zr^- Pr^)L + PlPlL^ (2.7) 
where the scalar L is defined as L = L ■ z in terms of the the four- vector [14] 

As it stands, formula (2.7) is of poor practical interest. It is essential to observe [1] that Z commutes with 
■ P^. Let us transform (2.7) in order to render this property manifest. First we split z as the sum of z^ 
and Zt in Z, hence 

Z = {zl + zl)P^ - {z^ ■ Pf - {z^ ■ Pf - 2{z^ ■ P){z^ ■ P) (2.9) 
Develop (2.7) and perform elementary manipulations using (2.9). We get 

Z = Z + 2Pl{z^ . P)L + 2Pl{z, ■ P)L - 2^{z, .pf + ^{z. P,f 

Using (2.9) again we obtain 

Z = Z^ 2{z^ ■ P){z, ■ P) + 2{z, ■ Pf - 2^{z, ■ Pf + ^{z, ■ Pf 



That is 



Z = Z + 2{z^ ■ P){z, ■ P) - {z, ■ Pf ^ 



Using (2.9) we notice cancellation of the terms proportional to {zt ■ P){zl ■ P) and we can write 

Z = zlP^ + zlP^ - {z^ ■ Pf - {z, . Pf{l + 



That is to say 

Z = r2 p2 _ . p^2 j_ p2(^2 _ (fi 

P? 

It is convenient to define the projector "orthogonal" to P^, say 



Z = zlP' - • Pf + P\zl - ^^^1^) (2.10) 



Q^ = <5|-:^^ (2.11) 

because we can write 

- ^^Sf^ = ^2.12) 
and we easily check that {Q.z)"' commutes with {y^ ■ Pl)- So we finally have 

Z = zlP^ - {z^ ■ Pf + {nzr.fP'^ (2.13) 
which justifies the claim that Z commutes with Hl ■ Pl Here we notice that VLZt = Zt and finally obtain 

Z = (nzfP^ - {z^ ■ pf (2.14) 

This simplification of (2.7) renders the ansatz more tractable. 

Remark: Formula (2.14) was derived without specifying the respective dimensions of the longitudinal 
and transverse spaces. It is valid also in the case considered in ref.[9], but in that case Q.Z reduces to Zt, 
which makes ^Iz^^ to vanish. 

It is convenient to replace the basic wave equations by their sum and difference, setting 

After transformation (2.4), and in view of (2.6) we need the expression of K[ + K'^. Equations (3.36) of 
Ref.[l] yield in the present notation 

K[+K'^ = Kt+K2- 2T{L ■ y) + T\L ■ L) (2.15) 

where 

T = A'l - K2 -yL-PL (2.16) 
and the four-vector L is given by (2.8). Wc notice that 

^•^=(i^' ^•^=^ ^'-''^ 

For neutral systems, a further transformation inspired by the work of Grotch and Hcgstrom [4] permits 
to get rid of the Q variables. We introduce a new wave function \E'" = (exp iF)^' with the help of the unitary 
transformation generated by 

T=^{z.F.Q) (2.18) 

We set 

= exp(ir) O exp(-ir) O" = {O'f yO (2.19) 
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By (2.19), pscudomonicntuni is transformed to P". Since i/l ■ Pl and P" commute with F, they are not 
affected by transformation (2.19), thus finally 

H'^ -H^ = K'l -K'i = y^-P^ (2.20) 

Taking (2.20) into account, consider sum and difference of the wave equations. 

{H'{ + H'^)^" = 11^" (2.21) 

■ = J/*" (2.22) 

Transformation (2.19) ensures that H'l + H!^ doesnot depend on the external variable Q. Moreover H[' + H2 
commutes with yL- Pl', therefore it should not involve the operator ■ Pl conjugate to it. We shall explicitly 
check this point in Section 4. 

From now on we demand that ^P" he eigenstate of pseudomomentum with a timelike four-vector k" as 
eigenvalue. Combining this requirement with (1.22) we obtain 

= exp(ifc • Q) expii^%r^) <l> (2-23) 

\kL\ 

where ^ depends on z, but only through its projection orthogonal to fci,, and additionally depends on k and 
on v as parameters. In other words (p — 4>{^i k, wz) where we define xu as the tensor which projects any 
vector on the 3-plane orthogonal to fci,. The linear space of such functions includes as a subspace the Hilbert 
space 

L^{kL)^L^{^^,d'^wz) (2.24) 

characterized by convergence of the triple integral / (j)*(l)d^wz which simply reads / (j)*(j)d'^z in any frame 
adapted to k^. We use the subscript k for typographical simpHcity; it refers to the vector k. In fact, the 
Hilbert space defined above depends on k only through the direction of kL (it doesnot depend on kf). 

Remark: There is no explicit contribution of external field to equations (2.22) (2.23). In contradistinc- 
tion, equation (2.21) remains sensitive to the presence of external field. 



3. EXPLICIT FORMULAS 

Our goal is to discuss the system (2.21)(2.22). In order to have an explicit expression for (2.21), let us 
first calculate H'^ + H!^. According to the notation set in (2.19) and to the definition of ^" we can write 

H'{ + H'i = {K[ + K'^f + 2V'^ = K'{ + K'i + 2V" (3.1) 



We have to transform {K[ +K2) and V as indicated in (2.19). But V is given by (2.6). According to (3.1), 

we need to transform formulas (2.15) and (2.6). Notice that the transformation generated by T always can 
be carried out explicitly. It leaves qi , 52 unchanged whereas 

Pi=Pi-|^-92 pl=P2 + ^F-qi (3.2) 

Hence 

P* = P+^F-z (3.3) 

y^ = y-lF-Q (3.4) 

Notice that F ■ z has only transverse components and they depend on the transverse variables Zt only. 
Therefore we can write 

P| = Pl (3.5) 



Eq (3.3) also implies 



pi^ = P^ + eP.F-z+^{F-zf (3.6) 
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{Plf = P^ + eiP-F.z) + ^{F. zf (3.7) 



These formulas will be useful later on. 

Straightforward calculation made in ref [1] yields 



Ki = \{p,-^-z.Ff (3.6 



Obviously = and was given in (3.3) above. Explicit development gives 

2 

2^1 =pI- |pia^ai^"" - |^ai^">la + ^{z ■ Pf (3.9) 

2 

2i^l =pI- lp2aZaF^" - |^<.F->2a + ^{z ■ Pf (3.10) 

We observe that all terms quadratic in the charge disappear from the difference; we get 

2{Ki - K2)" =vl-vl- eiVaZa + Zaya)F''" 

As [z„, Ha] is symmetric as a tensor, the order of y, z is immaterial when multiplied by the skew-symmetric 
F. Recall that Pi— P2 = V ■ P, hence the useful formula 

{Ki- K2Y =y P-2e z- F -y (3.11) 

But our goal was to transform eq (2.15). We first calculate {Ki + K^)"^. 
Notice that P ■ {z ■ F) = {z ■ F) ■ P = z ■ F ■ P, hence 

2{K^+K2f =pl+pl-ez-F-P+^^{z-Ff 

But pI+pI = P'^/2 + thus 

[Kr + K^f = ^ - ^-z. F . P +^^{z- Ff (3.12) 

Then we have to transform T, L ■ y and L ■ L. Take (2.16) into account and remember (3.11). (we know 
that y^ ■ Pl and -Pc are not affected by tt). We get 

T"^ = y^-P^-2ez-F -y (3.13) 

And finally after a glance at (2.17) we observe that L ■ y and L ■ L are unchanged in the transformation 

generated by T. 

Now we apply transformation (2.19) to (2.15), taking (3.12)(3.13) into account. It gives 

K'l + K'^ = 

^+y^-^z-F-P+^{z-F)^- 2{y^ ■P^-2ez-F-y)L-y+ {y^ ■ P^ - 2ez ■ F ■ yf L ■ L (3.14) 

We know that 2V" must be added to this expression in order to obtain H'{ + H2. But in (2.18) F^" is purely 
transversal, therefore (y^, • P^)' = y ■ P^. We have by (2.6) 

V" = f{Z\P^\y,-P,) (3.15) 
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where P'^ is as in (3.6) and we must compute from (2.13) by help of (3.3)(3.5). We make the convention 

that = (Z)" (and not the reverse). 

Then we apply the transformation (2.19) to eq. (2.13). Inspection of (3.5) shows that {^z^)^ is not 
affected by the transformation. We notice that Zt ■ = Zj. ■ P because, F being purely transverse, Zt - F ■ z 
identically vanishes. Thus 

= P^''{^lzf -{z^-Pf (3.16) 
Now, eqs (3.1)(3.14)(3.15) supplemented with (3.6)and (3.16) furnish the complete expression of H'{ + H!^. 

4. THREE-DIMENSIONAL REDUCTION 

4.1 The reduced M^ave equation for isolated systems 

Our presentation of the three-dimensional reduction will be more transparent if we first remind how it 
is usually carried out (and how it leads to a spectrum for P^) when the system is free of external forces. 

In this subsection wc drop the superscripts (0) referring to an isolated system. In the absence of external 
field, the squared-mass operators can be written 2Ha = Pa + 2^. A little of elementary algebra provides 

Hi + H2 = ^ + y^ + 2V (4.1) 

H,-H2 = yP (4.2) 

and the wave equations take on the form 

p2 

(— + 2/' + 2y)* = (4.3) 

[y . P)^r = (4.4) 

Introducing as defined in Section 2 we get 



p2 (v ■ P)"^ 

i?i+if2 = — + ^^^+r + 2F (4.5) 



and the wave equation (4.3) now reads 



Introduce the operator [6, 15] 
equation (4.6) takes on the form 



p2 (-„ . p^2 

(^ + ^^p2^+y' + 2m = M* (4.6) 



N = y^ + 2V (4.7) 



P^+^)* = '"* (4.8) 

It is noteworthy that —N is intimately related with the energy of relative motion (see the nonrelativistic 
limit A <C /i). Taking (4.1)(4.2) into account it is easy to check the identity 

_N=^+iJ[l^_^H,+H2) (4.9) 

Let us assume that 'J' is also eigenstate of the linear momentum, that is P^vp = /;-"^', where k is timelike. 

Since P^ and y ■ P are simultaneously diagonalized, so are the quantities P • pi and P ■ P2- It is usually 
required that the eigenvalues of these operators are both positive, which amounts to say 

Ifc^ > W\ (4.10) 
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Justification for this condition is easy when the mutual interaction is = f{Z,P^). In the imdcrlying 
classical theory [15] the equations of motion imply that P ■ z = P ■ p\ t\ — P ■p2 T2 + const, where n, T2 are 
the evolution parameters (generalizations of the proper times). On the equal-time surface (S) ~ P - 2: = we 
get the relation P-pi ti ~ P-p2 T2 + const. between the parameters. It is clear that ti must increase together 
with T2, which implies that classically P ■ p\ and P ■ p2 have the same sign (necessarily positive, since their 
sum is P^). Equation (4.10) is just the quantum conterpart of this condition. Extension of this argument 
to more general interactions is an open question, but other considerations motivated by the Bethe-Salpeter 
equation and propagator theory [16] strongly support condition (4.10). 

We can obviously write 

* = exp(ifc • Q)'4){z) (4.11) 
Using also (4.4) we cast eq. (4.8) into the form 

4 + ^+^)* = M* (4.12) 

Let us set 




(4.13) 



(4.12) now reads 

{N + \)^ = Q (4.14) 

which suggests to consider (4.14) as an eigenvalue problem where A is the eigenvalue for the operator —N. 
Caution is needed however, in the cases where V actually depends on P^ {energy- dependent case); in this 
case may be substituted to P^ into the potential and further eliminated from (4.14) with help of (4.13) 
(as we shall see there is only one way to solve (4.13) for fc^). But one is left with an extra dependence on A in 
equation (4.14). This complication, first pointed out by Rizov, Sazdjian and Todorov [12], will be discussed 
below, in the context of the reduced problem. 

Indeed equation (4.4) implies that the form (4.11) of the wave fimction can be further reduced. Recall 
that Ua = —id/dz"". In a frame adapted to fc, we can write y-k = j/o^" thus (4.4) entails ip = exp{ii'z^ /k°) 4> 
where doesnot depend on z*^. In fact 4> belongs to a linear space which depends on the direction of the 
vector fc", but doesnot depend on the number k — ^/k'^ka- Imposing square integrability defines a Hilbert 
space, say L^{k°'). Now eq. (4.8) or (4.12) becomes an equation to be solved in L'^{k°') because after 
substitution of k for P and v for {y ■ P) in N we obtain an operator acting in L^(fc"), say {N)ij,k- Insofar as 
V depends on P^, the operator {N)^^k depends on the scalar fc^. We can solve (4.13) for k and substitute 
the result in {N)v,k- In other words, we insert 

e = 2{\ + n)+2^{\ + -v-^ (4.15) 

into {N)^^k- 

Notice that an alternative solution of (4.13) corresponding to the minus sign in front of the radical is 
discarded because it would violate the condition ifc^ > \v\. In contradistinction, provided that A + ^ > 
we are sure that this condition is satisfied by solution (4.15). As a result the correspondance between A and 
fc^ remains one to one for unequal masses. 

But, for energy-dependent interactions, we end up with a reduced wave equation which fails to be an 
eigenvalue equation in the usual sense, being nonlinear in the eigenvalue A [17]. This complication, pointed 
out in ref.[12] stems from the possible dependence of the mutual interaction term V on the (squared) total 
energy P^. 

Fortunately, a mathematical theory exists for generalized eigenvalue equations which are nonlinear in 
the eigenvalue [18]. But this theory involves technical sophistications resorting to the concept of " associated 
vectors" . Its main departure from the standard theory is that eigenvectors corresponding to different gen- 
eralized eigenvalues may fail to be orthogonal. This drawback (not characteristic of the two-body problem, 
see [19]) led the authors of ref. [12] to advocate a re-definition of the scalar product of eigenvectors. In this 
paper we shall focus on situations where this redefinition is unnecessary. 



9 



No such problem arises when V docs not depend on the energy. Any potential of the form V = 
f{Z/P^, y ■ P) (in particular the harmonic oscillator V = const.z^) leads, in absence of external field, to an 
ordinary eigenvalue equation in L'^{k°') = L'^{'R^ ,(P z±), where z± is orthogonal to k. 

This class of potential encompasses a lot of phenomenological models. In contrast, realistic potentials 
motivated by QED pertain to the other type. 

4.2 The reduced equation in the presence of external field. 

Let us now return to the case where an external potential is applied to our system. The aim of this 
section is only to reduce the wave equations. The discussion of all matters concerning the spectrum are 

postponed to Section 5. 

Since C is the natural generalization of linear momentum in the presence of external field, and by 
analogy with (4.9), we shall define 

-N=^+^J^l^-iH,+H,) (4.16) 

Although formula (4.7) is no longer valid, the sum of the basic wave equations still reads (A'' + A)^' = 
with A defined as in (4.13). We assume that 'i' is eigenstate of the pseudomomentum, say C"'^ = k"^, 
for some constant k which is supposed to be strictly timelike. By analogy with the isolated system and 
also for pragmatic reasons, we shall confine ourselves to the sector k'^/2 > After transformation to the 
convenient representation we get 

p2 ( TTit _ fj>n2 

-N" = T + P2 - (^1 + ^2 ) (4.17) 

and our goal is to determine the spectrum of this quantity. Then one will pass from A to k"^ through (4.13). 

The calculations can be organized as follows: 

Whereas (2.22) fixes the dependence in the relative time, eq.(2.23) allows us to factorize out the "center- 
of-mass motion", and we are left with the reduced wave function which arises in eq. (2.23). Obviously 
(2.22) impUes that 

y^-k^0 = O (4.18) 

thus (j) depends on z only through its projection vjz. Imposing square integrability amounts to require that 
is in the Hilbert space L\ defined by (2.24). It is clear that generally depends on v and k as parmeters. 

In order to write down the reduced wave equation, we replace and • Pl respectively by their 
eigenvalues k" and v in iff + H2, and we divide by exponential factors. 

For any operator O it is convenient to use the following convention 

yL-PL=v, P=k (4-19) 
The subscript k refers to the vector k, which finally contributes by its longitudinal piece only. Notice that 
here a term like must be written as y^ = (fly)^ H ~"52~ — ■ ^® introduce the projector zu 

orthogonal to fci,, and use the identity -j = t5'(1 ~ e^) we obtain 

(^ + y\,k = ^ + iruyf + ^ = ^ + (-^)' + ^ - (4-20) 

which is to take into account when computing (AT" + K!^)^, !- from (3.14). 
Naturally {H'^ + H'^),,k = (AT" + A'^O^'.fe + 2{V").^k- Defining 

R{iy,k^,k^) = {K'; + K'^)^,k (4.21) 

W{iy,k^,k^) = {V'%,k (4.24) 
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we end up with the equation 

R4> + 2W 4) = ^i(t> (4.23) 

At this stage, let us observe that, in agreement with a remark made in Section 2, it is expected that neither 
R nor W involve z,^- P^^. This will be checked later on and will permit us to consider R and W as operators 
acting in the Hilbert space L/^^ki^) defined by (2.24). 

We remind that /x is just a parameter fixed from the outset. The question wether (4.23) can be con- 
sidered as a spectral problem, and for which eigenvalue, will be considered later on, with help of equations 
(4.13)(4.16). See eq. (4.33) below. 

The explicit expression of R comes from (3.14), with help of (4.21), and at first sight seems very 
involved. This apparent complication results from having used the necessary transformations (2.4)(2.19). 
Things become more tractable if we separate all terms involving F^^" from those which survive when the 
field is turned off, and so on. So let us perform such a separation in a systematic way. Be cautioned that 

(e)")(o) ^ (o(o))". 

Since K'{, K2 are no more than quadratic in the field strenght, we can write 

K'l + K'i = {K'l + K!^Y°^ + {K'l + K'if^^ + {K'{ + K'^f^^ (4.24) 

with the obvious convention that superscripts (1), (2) respectively refer to (homogeneous) linear and quadratic 
terms in the field strenght. Therefore R = + i?^^' + R'^'^\ In eq.(3.14) we have to replace P by A; and 
■ Pl by V, in order to compute R. Remembering (2.17) we start from (3.14) and compute K'{ + K2 (to be 
inserted into (3.1) and further simplified by the above substitution). We first consider the piece of it which 
survives in the no-field limit. In eq. (3.14) we collect zeroth order terms in the field and find 

{K'l + K'^r) = ^ + + {nyf + ■ p^Vll^l^y^llL (4.25) 

According to notation (4.19) we can write 

At this stage, one might be tempted to separate out the first two terms in the right-hand-side because they 
are constant. However the second one depends on k^^ which looses intrinsic meaning in the no-field limit. 

Therefore we prefer to apply the identity = -77(1 — e^) which yields 

kf k-' 

{K'l + = ^ + g - e^^ + i^yf + - 2.^ (4.27) 



Now setting 



we can write 



so we end up with 



j/t • fcy — 2j/ 2 '^^ 



{S)^,k = {wyY + y^ ■ ^^-^^^ 

(^f + ^2)S = ^ + ^ + (^k. (4.29) 



2 ,,2 



The field-depending terms in (3.14) provide 



R^^'> =4e{z-F-y)^-^z-F-k (4.31) 
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i?(2) = ^ (2 . + (^•-^^•^)' (4.32) 

Remember that F has transverse components only. Contractions involving F only depend on the transverse 
components; for instance F • fc is just a combination of the quantities fc". It is noteworthy that only the 
transverse components of z,y arise in i?'^) , whereas {S)i,^k depends on zuy and y^- As a whole, R 

depend only on zuz and zuy (recall j/t, Zt are pieces of zuy, zuz respectively). 

Let us again define A by formula (4.13). In view of (4.30)(4.31)(4.32), equation (4.23) may be written 

+ [{S),,k + R^^^ + -R*'^ + 2W](^ = (4.33) 

The square bracket in (4.33) is nothing but {—N")^^. 

Remark: R^^\ R^^^ R^'^^ do not involve the mutual interaction. In contradistinction W is defined through 
(4.24) and is model dependent in this sense that it crucially depends on the form of the function / which 
determines the mutual interaction. 

Owing to the abundance of terms depending on e we distinguish the "motional" case, where e ^ 0, 
from the case "at rest" characterized by the vanishing of e (or cquivalcntly of kr). Notice that R^'^^ doesnot 
depend on fc^ . It has the same form in motional case and in the case at rest [20] . 

Let us evaluate W. In view of (4.24) we have first to write down the expression for V" , say (3.15). It 
follows that 

W = f{{zXk, iP%k, (4.34) 

In this formula (P'')^ is given by (3.6) and by (3.16). Making the substitutions P — > and j/i, • Pl ^ f 
we obtain 

(^«).,fe = {P%k i^zf - (z. • kf (4.35) 

^P^\,^ = e + ek-F-z+^^{F-zf (4.36) 

It is clear that W docs not involve the operator z -k^. Formulas (4.35) (4.36) arc to be inserted into (4.34), 
then the explicit form of W will come out. This last expression, together with (4.30)(4.31)(4.32) provides 
the explicit form of the eigenvalue equation (4.33), where A is the eigenvalue of an operator —{N)^^k acting 
in 

Still we meet a complication: just like in the case of an isolated system, we can solve (4.13) for fc^ 
(condition (4.10) ensures that relation (4.13) can be uniquely inverted) and insert the result into {N")i,k- 

Similarly fc^ can be replaced by -. As a result eq. (4.33) bears an extra dependence on A. In general 

it is a nonconventional eigenvalue equation. 

Naturally a scalar product for the reduced wave function has to be explicitly defined. This question will 
be discussed in next section, in analogy with the line we have followed in the case of an isolated system. 

According to relation (4.13), a discrete spectrum for A would imply that the admissible values of fc^ are 
restricted to a discrete sequence. 

4.3. Discussion 

Finally the system (2.21)(2.22) has been reduced to the three-dimensional problem of solving (4.33). 
This formula is nonlinear in the field strenght and may be applied to strong fields. Let us analyze the 
various contributions it contains. We distinguish motional terms, depending on e or fcy- In fact we can write 
fcy = eAfci, where the second rank tensor A represent the boost from the direction of fci, to the direction of 
fcr (thus A • A = 5). 

Loosely speaking we could say that, in as much as the shape of W departs from the original form 
assumed by V^^\ the mutual interaction is "modified by the magnetic field". 

a) system at rest 

The particular case where pseudomomentum is purely longitudinal enjoys a particular simplicity. If we 
assume for a moment that k coincides with k^, it is possible to find a frame where k vanishes whereas the 
electromagnetic field is purely magnetic. In other words the pseudo rest frame is also a lab frame. We refer 
to this situation as the case at rest. 
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In this case, voz = z±, wy = y± and {S)^^k simply reduces to y'j^. Naturally L'^{ki^) coincides with 
L'^{k). We notice that the second term in the r.h.s. of (4.35) vanishes, hence reduces to z\. 

Thus, for energy independent interactions ( V^'^^ = g{Z/P^,y ■ P)), W assumes the form g{z^, v). In other 
words: At rest, the magnetic field doesnot modify the mutual interaction when it is not energy- dependent. 

All surviving terms in (4.33) can easily be identified as covariant generalizations of the usual terms 
present in the non-relativistic theory, except for a piece of i?^^' which depends on the relative angular 
momentum (formula (4.32)) and remains small for heavy bound states; its contribution however might be 
significant for lignt bound states (fc^ small) in a strong magnetic field. 

At first order in the field strenght, the relative motion admits no correction other than a term propor- 
tional to u (indeed F ■ k vanishes). For equal masses there is no departure from the motion of an isolated 
system. 

b) motional case 

When is nonzero, we reckognize the motional electric field contained in _F • fc. 

Energy-dependent interactions seem to be more sensitive to the external field. As can be read off from 
(4.36), the modification implied by (4.34) is nonlinear in F. This point may become important in strong 
fields. 

Even if the interaction is independent of the energy, the presence of the motional term z^. ■ k in (4.35) 
entails a non-trivial difference between (4.34) and (2.2). 

5. NORMALIZATION AND PERTURBATION THEORY. 
5.1. Normalization 

In connexion with the rise of a reduced wave function (p in (2.23), we are led to consider in general any 
function (or possibly a distribution) of z of the form (p{k, z) which depends on z only through the combination 

vjz (assuming k timelike), regarding the four- vector A: as a given parameter. The most straightforward 
normalization of 4> is given by a three- fold integral using the volume element d^{wz).i with the convention 
that wz = z in any frame adapted to /c^. In other words, for each fc, we separate fci, from k and consider 
the Hilbert space L^ik^) = L'^ {K^ , d^ {zu z)) . 

For a function of the above type, but not necessarily solution to (4.33), we have, with an obvious notation 

<(f>,(l>>k= j <f>* (t>d^{wz) (5.1) 

Here may additionally depend on v and kr as parameters. The label k refers to the vector k, but actually 
L'^{k^) depends on fc" only through the direction of k^- 

It is noteworthy that the normalization (5.1) of reduced wave functions is not only dictated by simplicity, 
but also consistent with the off-shell normalization of the full wave functions. 

The proof of this point will be given elsewhere; its precise meaning is as follows: 

Let us consider any fmiction "^"(Q, z) irrespective of its being on the mass shell or not. 

Assuming that (j){z°',k^, v) is related to through (2.23), it is possible to check that 

j d'^Qd'^z = j <(t),(l)>k d'^kdy 

provided that the Fourier transform of \E'" with respect to Q, say 5(fc,z), is a retarded function of k (it 
vanishes outside the region limited by the positive sheet of the light cone). 

Of course this nice property [21] doesnot prevent the complications associated with energy dependence. 

Even in the simple case where the mutual interaction (present in the isolated system) doesnot 
depend on P^, we are bound to realize that the spectral problem is not a standard one, due to the presence 
of the magnetic field. This can be seen as follows. First we notice that the occurence of {Z'^)u,k in W 
brings out a dependence on k'^ ,k'^. Second we observe an unescapable dependence on k^ and /c^ in formulas 
(4.28)(4.31)(4.32). Clearly, we cannot expect a conventional spectral problem. This peculiarity arises because 
going from ^ to we have abandoned the customary representation; transformation formulas (2.4)(2.19) 
are responsible for the nonlinear dependence on A. 
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5.2 Outline of a perturbative approach. 

In order to solve the reduced wave equation, perturbation theory is by no means straightforward. The 
abundant mathematical literature devoted to nonlinear eigenvalue equations [18] pays very little attention 
to spectrum perturbation. So the general case seems to be hopeless in the present state of the art [22]. 

Henceforth we limit ourselves to the simple case where F'-"-' is not energy dependent. Still we must cope 
with a nonconventional eigenvalue equation when magnetic field is present. But we can make a couple of 
important remarks. 

First we observe that all terms involving energy dependence vanish when both the external field and the 
motional parameter vanish. 

Second, we notice that e and k'^ are two independent parameters contributing to the determination of 
the vector fc". Inspection of (4.13) shows that e and A are mutually independent. 

These remarks open the possibility to treat simultaneously electromagnetic contributions and motional 
terms as a perturbation. This procedure requires that the motional parameter is not too large. 

In this approach the unperturbed equation is free of nonlinear dependence on A. It corresponds to 
"static states" of the isolated system, i.e. states with external coupling removed and additionally at rest 
with respect to the special lab frame; they are characterized by fc^ = 0, F — (sec Appendix). Only 
the perturbation, which combines motional terms and field contributions, is affected by energy dependence. 
It is fortunate that Rizov, Sazdjian and Todorov [12] have developed a perturbative scheme (tractable by 
physicists) for this situation. 

Another departure from conventional perturbation theory lies in the fact that the perturbation is non- 
linear in e and in the field strenght. 

Indeed, insofar as / in (2.6) can be developed as an analytic function of its arguments, a look at 
(4.35)(4.36) indicates that even if / is linear in Z (harmonic potential), V" and W are at least quadratic in 
the external field. 

Fortunately, the formalism set up in [12] is general enough to accomodate second (and higher) order 
perturbations. 

The first step is a weak-field-slow-motion approximation. The field strenght and e are taken into account 
only at first order. The perturbation scheme of ref. [12] gets simplified; it turns out that the correction to the 
eigenvalue is given by the usual formulas. In particular if the unperturbed eigenvalue is nondegenerate, the 
correction is still given by the expectation value of the perturbation term in the eigenstate of the unperturbed 
problem. We can organize calculations by setting F/ ji = eG (or alternatively F/k"^ = eG) where the tensor 
G is dimensionless. It follows that F ■ k and the motional electric field are O(e^). The first order treatment 
entails enormous simplifications. For instance the obligation to distinguish iz±)'^ from (tuz)^ is in general 
a serious computational complication, but these quantities coincide at first order in e. A glance at (4.35) 
shows that W = f{zj_, v) + O(e^). At first order, the only departure from the unperturbed equation consists 
in the i^-terms and vanishes for equal masses. To summarize: 

// the neutral system as a whole undergoes a slow motion in a weak field, the relative motion is affected 
by the magnetic field only by a contribution which vanishes for equal masses. 

Therefore it is necessary to go beyond first order in search for a nonzero correction. 

Example: 

Suppose that F^^^ = g{z^). This form of mutual interaction includes the relativistic harmonic oscillator, 
say = jZ/P"^, where 7 is a coupling constant. The only difference from the unperturbed equation comes 
from _R(^\ which is given by (4.31). 

Another situation of great simplicity consists in a system of equal masses at rest (e and u are strictly 
zero). The perturbation reduces to R^'^\ hence is linear in Corrections to the levels are easy to calculate, 
but significant only in strong fields. 
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6. CONCLUSION. 

The coupled Klein- Gordon equations describing a neutral system have been reduced to a three-dimensionall 
eigenvalue equation involving truly motional effects and recoil effects in a covariant fashion. Moreover the 
particular symmetry associated with a constant magnetic fic;ld in space-tim,6 is manifestly respected. 

Explicit formulas are written in a representation which permits to easily satisfy compatibility and also 
to eliminate relative time and the Q- variables. The surviving degrees of freedom are the same as in the 
nonrelativistic theory. 

In the reduction procedure it was essential to consider eigenstates of pseudomomentum. The square of 
this vector plays the role of an effective squared mass which can be, in principle, evaluated by solving the 

reduced equation. Indeed (in the sector we have considered, and especially for equal masses) the eigenvalue 
A is in one-to-one correspondance with k'^. Bound states are characterized by A in the discrete spectrum. 
In the most natural way, our approach generalizes to the magnetic case the usual treatment of isolated 

two-body systems according to "predictive mechnics" and "constraints theory" [6] [7] [16] [8]. It provides a 
clean theoretical basis for the study of neutral bound states in the relativistic regime. 

The equations of motion are nonlinear in the field strenght and offer a starting point for investigation 
of strong field effects. In principle, they encompass all kinematic possibilities of the system as a whole 
and permit a description of ultra-rclativistic situations. The Ansatz which allows for a three-dimensional 
reduction in the covariant framework automatically generates various terms in the wave equation. From a 
practical point of view, some of them play the role of a modification of the mutual interaction, caused by the 
magnetic field. This effect might have dramatical conseqiicnces in a system where the mutual interaction 
is energy-depending (for instance Coulombian interaction along the lines of the quasi-potential approach). 
But in this case we cannot go beyond qualitative estimations, since the perturbative scheme borrowed from 
ref. [12] is limited to the systems where the interaction is not energy dependent. This point may be 
a motivation for eventually undertaking a more general perturbation theory applicable also to eigenvalue 
equations involving an energy-dependent potential. 

When does not depend on P^, spectrum calculations can be handled in a perturbative approach 
where the motional parameter as well as the field strenght must necessarily be small. For the moment this 
method can be used for simple and idealized systems bounded by a phenomenological potential. Application 
to the relativistic "naive quark" model of hadron seems to indic;ate that the spectrum enjoys some kind of 
stability with respect to magnetic perturbation. Of course the polarizability of this model can be studied in 
our formalism. 

In the hope of obtaining a nonvanishing correction due to the presence of external field, we are obliged 
to go beyond first order. In a future work we plan to investigate with more details the case e = u = 
sketched in last section. 

More work is needed in search for ultra relativistic effects ( e ~ 1). A less ambitious program might be 
the computation of second order motional effects in the framework of the perturbative scheme considered 
here. But the complication of the calculations may be prohibitive. 

Naturally an extension to particles with spin is desirable. 

Let us finally mention that the contact with usual methods of quantum field theory could be improved 
by considering a Bethe-Salpeter equation taking the magnetic field into account from the start. This would 
mean to resume the work of Bijtebier and Broeckaert [23] in a way which respects the particular symmetry 
of constant magnetic field, i.e. treating all the possible lab frames on the same footing. 



APPENDIX 

The unperturbed equation 

The vanishing of kr (equivalently of e) makes L'^{ki^) to coincide with L'^{k) and wz with z±. As also F 
is zero, inspection of (4.34)(4.35)(4.36) shows that W reduces to f{{Z)^^k,k'^,i') = f{z'^k'^ - {z ■ k)'^,k'^,u) 
Finally we see that, for e — F — 0, (4.33) reduces to 

A0 + zlcj) + f{z^k^ - {z ■ kf, P, ly) (j3 = 

It is exactly the equation one would obtain for the isolated system after reduction, in the original represen- 
tation. 
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